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Abstract 

£jq[ The recent observations imply that there is an extra relativistic degree of freedom coined dark 

radiation. We argue that the QCD axion is a plausible candidate for the dark radiation, not 
only because of its extremely small mass, but also because in the supersymmetric extension of the 
Peccei-Quinn mechanism the saxion tends to dominate the Universe and decays into axions with a 
sizable branching fraction. We show that the Higgsino mixing parameter /x is bounded from above 
^■i when the axions produced at the saxion decays constitute the dark radiation: /j, < 300 GeV for a 



saxion lighter than 2mw, and fi less than the saxion mass otherwise. Interestingly, the Higgsino 
can be light enough to be within the reach of LHC and/or ILC even when the other superparticles 
are heavy with mass about 1 TeV or higher. We also estimate the abundance of axino produced 
by the decays of Higgsino and saxion. 
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I. INTRODUCTION 



The present Universe is dominated by the dark sector, i.e., dark matter and dark energy, 
although it is not yet known what they are made of. Therefore, it may not be so surprising 
if there is another dark component, which behaves like radiation. 

The presence of additional relativistic particles increases the expansion rate of the Uni- 
verse, which affects the cosmic microwave background (CMB) as well as the big bang nu- 
cleosynthesis (BBN) yield of light elements, especially 4 He. The amount of the relativistic 
particles is expressed in terms of the effective number of light fermion species, N e g, and it 
is given by N e g m 3.046 for the standard model. Therefore, if iV e ff > 3 is confirmed by 
observation, it would immediately call for new physics. 

Interestingly, there is accumulating evidence for the existence of additional relativistic 
degrees of freedom. The latest analysis using the CMB data (WMAP7 [l] and SPT 0) has 
given N e g = 3.86 ± 0.42 (la C.L.) j^J. Other recent analysis can be found in Refs. [l], Q, U- 
8j. The 4 He mass fraction Y p is sensitive to the expansion rate of the Universe during 
the BBN epoch 1 , although it has somewhat checkered history since it is very difficult to 
estimate systematic errors for deriving the primordial abundance from 4 He observations 16 ]. 
Nevertheless, it is interesting that an excess of Y p at the 2a level, Y p = 0.2565±0.0010 (stat)± 
0.0050 (syst), was reported in Ref. 17| . which can be understood in terms of the effective 
number of neutrinos, N eS = 3.68^7° (2a). 2 Interestingly, it was recently pointed out that 
the observed deuterium abundance D/H also favors the presence of extra radiation 19|, |20| : 
iV eff = 3.90 ± 0.44 (la) was derived from the CMB and D/H data Q- It is intriguing that 
the CMB data as well as the Helium and Deuterium abundance favor additional relativistic 
species, AN e g ~ 1, while they are sensitive to the expansion rate of the Universe at vastly 
different times. 

The extra radiation may be dark radiation composed of unknown particles. Then it is 
a puzzle why it is relativistic at the recombination epoch, why the abundance is given by 
AiVefj ~ 1, and why it has very weak interactions with the standard- model particles. 



1 Y p is also sensitive t o la rge lepton asymmetry, especially of the electron type, if any [9- 15 1. 

2 The authors of Ref. [18| estimated the primordial Helium abundance with an unrestricted Monte Carlo 
taking account of all systematic corrections and obtained Y p — 0.2561 ±0.0108 (68%CL), which is in broad 
agreement with the WMAP result. 
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In fact, there is a well-motivated particle with the desired properties, namely, the QCD 
axion. The axion appears in the Peccei-Quinn (PQ) mechanism 211 ] . one of the solutions 
to the strong CP problem [22|, in association with the spontaneous breakdown of the PQ 
symmetry. The axion remains extremely light: its mass is originated from the QCD anomaly 
and is in the range of 10~ 5 eV < m a < 1CT 3 eV for the PQ breaking scale in the cosmological 
window, F a = 10 10 GeV ~ 10 12 GeV. Furthermore, in a supersymmetric (SUSY) framework, 
the saxion tends to dominate the energy density of the Universe, and it decays mainly into 
a pair of axions. Such non-thermally produced axions naturally remain relativistic until 
present. Thus, the axion is a plausible candidate for the dark radiation. 3 

In this paper we consider the QCD axion as a candidate for the dark radiation, and show 
that the branching fraction of the saxion into axions naturally falls in the right range, if 
the axion multiplet is coupled to the Higgs superfields as in the Dine-Fischler-Srednicki- 
Zhitnitsky (DFSZ) axion model 3l|, |32|. Most important, we find that the /i-parameter is 
bounded from above, namely fi < 300 GeV for the saxion lighter than 2mw, and \i less than 
the saxion mass otherwise, independent of the PQ scale F a . As we shall see later, if the 
radiative corrections play an important role in the stabilization of the saxion, its mass can be 
naturally smaller than the soft masses for the SUSY standard model (SSM) particles. This 
implies that, even if the other SUSY particles are so heavy that they are above the reach of 
the Large Hadron Collider (LHC), the Higgsino should remain within the reach of the LHC 
and/or future collider experiments such as International Linear Collider (ILC). This will 
be of great importance especially if the mass of the SM-like Higgs boson is indeed around 
124 - 126 GeV as suggested by the recent data from ATLAS and CMS at the LHC |33|. 

We also estimate the axino production from the decay of both Higgsino and saxion. In 
our scenario, the dark matter is made of the axion in the form of non-relativistic coherent 
field oscillations and the axino lightest SUSY particle (LSP), while the axion produced by 
the saxion decay accounts for the dark radiation. Thus, the axion and its superpartners play 
an important role to account for the dark matter and dark radiation. 



3 The abundance of relativistic axions produced by flaton decays was studied in detail in Ref. [23j . The 
late-time increase of N e g by decaying particles (e.g. saxion to two axions, and gravitino to axion and 



axino) was studied in Ref. 24 L The light gravitinos produced by inflaton decay can account for dark 
radiation in analogous to Ref. |25[. See also recent Refs. 2614291]. The possibility that the X particle was 
in thermal equilibrium was studied in Ref. [30 1 . 
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II. PQ EXTENSION OF THE MSSM 



A. Saxion Properties 

Let us begin by examining the properties of the saxion in a simple PQ extension of the 
minimal SSM (MSSM) where the PQ sector consists of the axion superfield S and pairs 
of PQ messengers ^ + ^ forming 5 + 5 representation of SU(5): 

W PQ = yySt>y. (1) 

Since the saxion is a flat direction in the supersymmetric limit, its properties are determined 
by how SUSY breaking is transmitted to the PQ sector. It is natural to expect that the PQ 
messengers, which are charged under the SM gauge groups, feel SUSY breaking in the same 
way as the MSSM superfields do. 

The PQ messengers radiatively generate a potential for the saxion after SUSY breaking 



341 ]. as well as an effective coupling of S to the gluon supermultiplet that implements the 



PQ mechanism. To see this, one can integrate out \1/ + \l/ under a large background value of 



5*, which results in 



£ eff = J d A 9Z s (Q = y*\S\)\S\ 2 . (2) 

There then arises a potential, V Ta .d = m2 s(Q — y^lSDlS] 2 . Here m| is the soft scalar mass 
squared of S, and Q is the renormalization scale. Hence, if m| is positive at a high scale 
and the messenger Yukawa coupling is large enough to drive it negative through radiative 
correction, Kad develops a minimum along \S\ around the scale where m| changes its sign. 
At the minimum, the Kahler potential fl5]) gives masses to the saxion and axino 



m 



a 



^2 Yl yl( m l + m l + A l) ~ 5 ^2 ff ' 

— ^y 9 A v ~ -^-Av, (3) 



where m 2 is a soft scalar mass squared, and A^ is the trilinear soft parameter associated 
with y^. The saxion acquires a mass suppressed by a factor e ~ A/5iV^/87r 2 y^ relative to 
other superparticle masses m so ft, and the axino has a small mass m„ ~ e 2 Aq,. From the 
effective Kahler potential, one can also find the saxion couplings 

u /,„., _ 1 



V2c a F a 



K^aJ^a + ^(A g ra 5 aa + h.c.H , (4) 



where F a = c a 1 (|S'|) is the axion decay constant. The model-dependent parameter c a is 
generally of order unity, and the saxion coupling to the axino 



(5) 

Q=c a F a 



dlnQ 

generically has a value in the range between 10 _1 and 10~ 2 . 

On the other hand, depending on the mediation mechanism of SUSY breaking, V ra d alone 
may not be able to stabilize S. This happens, for instance, in gauge mediation, where 
only those carrying SM gauge charges acquire soft SUSY breaking masses at the messenger 
scale M mess . An interesting way to give a soft mass to the gauge singlet S at M mess is to 
consider mixing between \l/ (or and the messenger fields transferring SUSY breaking 
35|. Then, there is an additional contribution to soft masses for the PQ sector fields, while 
the MSSM sector feels SUSY breaking only through the gauge mediation. The additional 
contribution makes V va & develop a minimum in the same way as discussed above. Turning 
off the mixing, the saxion potential runs away to infinity, but can be lifted by supergravity 
effects AV = ^m^ 2 \S\ 2 where m 3 / 2 is the gravitino mass. For a positive £ of order unity and 
m 3 / 2 not so small compared to m so ft, S is stabilized below M mess by V iax x+AV. The situation 
is different when m so f t 3> m.3/2, for which the supergravity effects become important at \S\ 
larger than M mess . 4 

We close this subsection by stressing that the potential generated by PQ messenger loops 
is naturally expected to play an important role in stabilizing the saxion in models where the 
PQ symmetry is spontaneously broken mainly by a single PQ field. Then, the saxion and 
axino are generally lighter than the MSSM superparticles. 



B. Saxion Cosmology 

Let us move to the cosmology of the saxion. The saxion has a very flat potential lifted 
only by the SUSY breaking effect, and so, the saxion potential may be significantly mod- 
ified during inflation, because SUSY is largely broken by the inflaton potential energy. In 

4 In this case, a minimum appears at a scale above M mess , where + W generate a potential for the saxion 
at the three-loop level [36j. This results in that the saxion has a mass of the order of TO3/2, and the 
axino acquires a tiny mass ma = (ei + £2)^,3/2 "C "1,3/2: where t\ ~ TO3/2/(87r 2 m so ft) comes from the PQ 
messenger loops while £2 ~ F% l^-pi i s due to the supergravity contribution. The coupling Aa has a value 
about (ei/871- 2 + e 2 )/(ei + e 2 ), and F a ~ M mcss TO so f t /(\/87r 2 "m3/ 2 ). 
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particular, the so-called Hubble-induced mass term generically deviates the saxion from the 
low-energy potential minimum. 5 

For instance, the saxion may be stabilized at a large field value close to the Planck scale 
in the presence of the negative Hubble-induced mass term. In this case, the saxion starts to 
oscillate with a large initial amplitude when the Hubble parameter becomes comparable to 
the saxion mass m a . Then the saxion would dominate the energy density of the Universe, 
if the reheating temperature of the inflaton is higher than the saxion decay temperature. 

Alternatively, the PQ symmetry may be restored during and/or after inflation. In particu- 
lar, since the PQ messengers generate a thermal potential for the saxion, thermal ~ V%T 2 \S\ 2 
at \S\ <C T, the saxion could be thermally trap ped at the origin if it sits around the origin 
after inflation. In this case thermal inflation [37J takes place when the saxion potential 
energy at the origin (~ m^F 2 ) dominates the Universe. 

Thus, it is plausible that the saxion dominates the energy density of the Universe. In the 
rest of the paper, we assume that this is indeed the case, and discuss how the saxion decay 
proceeds. 



C. PQ Solution to the [i Problem 

If S has no other interactions with the MSSM sector than the loop-induced coupling to 
the gluon supermultiplet, the saxion dominantly decays into axions with a branching ratio 
B a ~l: 

r> ^ c ¥0,0, _ 1 1 ml 

a r CT r a 6^(c a F a r [) 

where To- is the total decay width of the saxion. This is obviously problematic in a scenario 
where the Universe experiences a saxion-dominated epoch. 

A natural way to suppress B a is to introduce a coupling of S to the Higgs doublets so 
that a Higgs \x term is dynamically generated after PQ symmetry breaking. Among various 
types of couplings, we take the Kim-Nilles superpotential term 



W = \-^H u H d} (7) 
Mpi 



This is not the case if there is a certain (approximate) symmetry among the PQ fields. 
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which generates an effective ji term 



/i = A 




(8) 



Mpi 



around the weak scale for A < 1 and 



F a = 10 10 12 GeV, and thus naturally explains the 



smallness of /i. 6 The above superpotential term also induces couplings of the saxion to the 
MSSM particles. These couplings are proportional to fi/F a , and open up the possibility to 
obtain AiV e ff ~ 1 from the axions non-thermally produced by saxion decays. 

To evaluate B a , one needs to know interactions of the saxion with the MSSM particles. 
For simplicity, we consider the decoupling limit where the effective Higgs sector below m so f t 
consists only of a SM-like Higgs doublet and other heavy Higgs bosons decouple from the 
theory. Then, the saxion has the interactions 



both of which are directly induced from ([7]). Here h is the CP-even neutral Higgs boson, 
is the mass of the CP-odd neutral Higgs boson, and B is the Higgs mixing soft parameter. 
The coefficient C a = d\n\n\/dln\S\ = 2 reflects the form of the coupling between S and 
the Higgs doublets. In addition, the saxion interacts with other MSSM particles as well 
because there arises mixing between the saxion (axino) and neutral Higgs (Higgsinos) after 
electroweak symmetry breaking. These couplings can be read off from the MSSM Lagrangian 
by taking the substitution 



where (\H®\) = v sin/3. It is straightforward to derive the partial decay widths of the saxion 
into SM particles and Higgsinos, which are presented in the appendix. The important decay 
channels are those into hh, WW, ZZ, bb, ti, and into a Higgsino pair, depending on m a . 
We note that B a is determined essentially by fi and m a , but insensitive to F a . The saxion 
interacts with the SM particles more strongly for a larger /x, making B a smaller. 

6 To solve the fi problem, one may instead consider the Giudice-Masiero mechanism [3^ implemented by a 
Kahler potential term S* H u Hd/ S . In this case, however, the couplings of the saxion to the MSSM sector 
arise at the loop level. 






(10) 
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III. AXION DARK RADIATION 



In this section, we examine the condition for the axion dark radiation produced at saxion 
decays to yield AiV cff ~ 1 in the presence of the superpotential term (j7|) that is responsible 
for generating a p term. Then, we also examine if the axino can constitute the dark matter 
of the Universe with the correct amount of a relic density. 

Our analysis does not assume any specific mechanism of the saxion stabilization, but the 
radiative potential induced by the PQ messengers may play an important role in fixing the 
vacuum expectation value of S. Keeping this in mind, we shall consider the case where 
the saxion and axino have masses lighter than m so f t . It should be noted that the Higgsino 
mixing parameter p, on which B a crucially depends, does not need to be related to the scale 
?7tsoft because it is a supersymmetric coupling. 



A. Axion Dark Radiation from Saxion Decays 

The axions produced at the saxion decays have never been in thermal equilibrium for 
the saxion decay temperature much lower than F^/Mn J^, which is the case we shall deal 
with. The energy density of non-thermalized axions and thermalized radiation evolves as 

p a (t) oc a~~ 4 (t), 
Psu(t) oc g-y 3 (t)a- A {t), (11) 

where a(t) is the scale factor, and p a : psm — B a : (1 — B a ) at the saxion decay time t = t a . 
Because non-thermally produced axions contribute to N e g with 



AN eS = El 



Pv 



PSm(*) Pait) 

X 



v decouple 



v decouple p v (t) v decouple PSM (t) 

.H_^(i|i) ,/3 , (12) 
7 1 - B a \g*{T a ) J 

ANes = 1 is achieved if the saxion couplings to the MSSM sector arising from the superpo- 
tential ([7]) are strong enough to give B a between about 0.2 and 0.3. Here p u is the energy 
density of a single species of relativistic neutrino, and the effective number of relativistic 
degrees of freedom varies from ~ 60 to g* = 113.75 for 0.2 GeV < T < m so f t , through 
which a mild dependence on T a comes in. The saxion decay temperature T a is defined at 
the time t = t a when the energy density of radiation p a + p SM becomes equal to that of the 
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FIG. 1: The constant contours of AiV e g for \B\/mA = 0.6 and = 125 GeV in the (m a ,fj,) 
plane. The black lines represent the contours of AN e g = 0.5, 1, 1.5 from the above, respec- 
tively. In the shadowed region, 0.4 < A7V e ff < 2. We also plot constant contours of the quantity, 
(c a -Fa/10 11 GeV) x Tfj — 0.1, 0.5, 1, 3, 6, 10 GeV from the below, respectively, in red lines. 



saxion, which implies Psm(^o-) = (1 — B a )p a (t a ), where we have used that the ratio between 
energy densities of axion and thermalized radiation is p a /psM = B a / (1 — B a ) at t = t a . One 
can then estimate, 

1/4 



--mm ^ 



„ fB 2 J(l-B a )Y 1/4 fgJT a )Y 1/4 / m a \3/2 / Ca F a \ _1 . s 
. 3 .5Gev(^L_^) (^) (3^) ,(13) 

with t a ~ l/r CT , by taking the approximation that the scale factor a(t) of the Universe is 
determined mainly by the saxion energy density at t < t a and by the energy density of the 
radiation at t > t a . 

In fig. [H we show constant contours of AN e s in the (jn a ,fi) plane for \B\/mj\, = 0.6 
and rrih = 125 GeV. Here we take g„ = 80, which depends on T a but modifies AN& only 
slightly for 0.2GeV < T G < m so f t . It is interesting to observe that AN e s ~ 1 is obtained at 
p of a few hundred GeV for m a less than 1 TeV. A large \x renders the saxion couplings to 
the SM particles strong, suppressing B a . We find that non-thermally produced axions yield 
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AiVeff = 1 at /i similar to the saxion mass: 



fi ~ 150GeV x 



1 



\B 



-1/2 



GoOGev) ' (M) 



for 2my/ < m a < rn so f t , in which case the saxion decays mainly through a — > WW and 
through a — > ZZ, hh, ti if kinematically allowed. Here mw is the W boson mass. On the 
other hand, for m a < 2mw, where the dominant decay channel is o — > bb, a fi less than 
about 300 GeV is necessary to have AA e ff = 1. In the figure, we also draw constant contours 
of the quantity 

ira) T - (15) 

from which the value of T a can be read off for a given value of F a . As we will see soon, the 
dark matter abundance puts a stringent constraint on T a and also on the properties of the 
axino. 

The saxion generically acquires a mass smaller than m so ft by one order of magnitude from 
the potential radiatively generated by the PQ messengers. Meanwhile, the recent data from 
;he LHC Higgs search suggests that a SM-like Higgs boson may have mass around 125 GeV 



331 ] . To explain this within the MSSM, we need to push the stop mass to about 10 TeV 
or higher, or to invoke large stop mixing. For m a < 0.1m so ft and ITeV < m so ft < lOTeV, 
the Higgsino should have a mass of a few hundred GeV when the axion dark radiation gives 
AiV e ff ~ 1- Thus, there is a chance to detect SUSY at multi-TeV hadron colliders even when 
the other MSSM superparticles are heavier than 1 TeV. 

B. Axino Dark Matter 

Under the i?-parity conservation, the axino is a natural candidate for the dark matter 
in PQ extensions of the MSSM. The main processes of axino production for fi < m so f t are 
(i) the decay of Higgsinos in thermal bath, and (ii) the decay of the saxion, which crucially 
depends on the properties of the axino. 

Let us first examine the thermal process (i), which is mediated by the interactions 

with Ca = d In /i / d In S = 2, where the first coupling comes from the superpotential term 



-h (h° u cos (3 + H° d sin p) a + -jL^Z^ fe sin (3 - H° d cos ) a^fi, (16) 



S 2 H u Hfi, while the other is a consequence of the axino-Higgsino mixing. The produced 

10 



axino number density is highly sensitive to T CT , 7 and is numerically approximated by ria/s oc 
e -o.63^/T CT f Qr T/ < To- < 0.2/i where Tf ~ /i/20 is the freeze-out temperature of the Higgsino. 
The derivation can be found in the appendix. We find that n„/s < 10 -10 requires T a < Tf 
for fi ~ m a and -B a ~ 0.1. For T a <Tf, the axino energy density is approximated by 

/ dt a3 (*) r £? n £P (*) + n H« (*/) 



Pa 
S 



ma 

thermal S 



(17) 



where r& is the decay rate for H® — >■ h(Z)+a, n&o is the number density, and t/ is the time 
when the Higgsino freezes out of thermal equilibrium. We have taken into account that there 
are two neutral Higgsinos H® 2 , which are almost degenerate in mass, and neglected mixing 
between Higgsinos and gauginos. Meanwhile, because m a < 2/i is required for AN e g < 2 
unless B is much smaller than m so f t , the saxion decay into a Higgsino pair is kinematically 
forbidden. Even if this mode is open, the annihilation among the Higgsinos produced at 



the saxion decays would occur effectively 45), and consequently the axino relic abundance 



produced from the Higgsino decays can be smaller than the observed dark matter abundance 
for rria less than T a . 

On the other hand, the non-thermal process (ii) crucially depends on the properties of 
the axino. This process gives rise to 



Pa 
S 



2r CT _ > aa ?7t a 
non— th O^S Tfl a Jq 



dt a 3 (t)p a (t) 



3e r ^ T n m~„ T 



2 1 - B a rrirr T r 
2.4 x l(T 10 GeV ( ) ( v ^ ~ aj ) ( - ) ( ) . (18) 



Aa V / B a /(1 - B a )\ ( ma/m a \ 3 ( T a 



0.01 J V 0.3 J V 0.01 J \lGeV, 
Thus, in order not to overdose the Universe, the axino should have a small coupling to the 
saxion and/or small mass compared to m a . The non-thermally produced axinos can yield 
VLah 2 ~ 0.1, for instance, if m- a is less than a few GeV and \- a < 10~ 2 for F a = 10 10 ~ 12 GeV. 
Such axino properties are indeed expected when S is stabilized by the potential generated 
from the PQ messenger loops. 



7 As was noticed in Ref. 4l|, in the case that the Universe does not experience the saxion domination, 
thermally produced axinos by Higgsino decays would overdose the Universe unless the reheating temper- 
ature is much lower than the weak scale or the axino has a tiny mass less than 0(100) keV for F a ~ 10 11 
GeV. See also Refs. 



42 



44j for thermal production of the axino. 
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The axino abundance, (flTl) and ([TBI) , can be made consistent with the observed value of 
the dark matter abundance by taking appropriate values of m„, X a and F a , for fi and m a 
leading the axion dark radiation to give AN e s ~ 1. The energy density of the dark matter 
receives contribution also from the axion due to the vacuum misalignment 

F a \ 1,18 

lO^^evJ ' (19) 
where \9 a \ < 71 is the initial misalignment angle. If the axino relic abundance is too small, 
which would be the case for F a around 10 12 GeV or higher, the dark matter of the Universe 
can be explained by the axion from the misalignment. 

One might consider other cases where the axino is heavier than the Higgsino or some other 
MSSM superparticle. Then, for \i of a few hundred GeV as required to have AA^ ~ 1, one 
needs either m a < 2m a or an extremely small A„ in order to avoid overproduction of the 
dark matter. In the case where the gravitino is the lightest superparticle, which is possible 
in gauge mediation, a small gravitino mass m 3 / 2 <C m„ or a tiny A„ would be necessary since 
gravitinos produced at the axino decays behave like a hot dark matter with a free-streaming 
length much larger than 10 Mpc. 

Finally, we mention the detection potential of SUSY at collider experiments. The charged 
Higgsino generally obtains a mass slig htly heavier than the mass of when bino and wino 



masses have the same phase 46|, |47| . Assuming that it is the lightest one of the MSSM 



superparticles, H® decays into h(Z) + a with 

r .„ « ^l^l L_ ( H V ( ^ ^ ~ 2 m 

H ° 16vr(c a F a ) 2 " 314cm V200GeW Vl0 11 GeVy' ' 1 ' 

for fi > rrih + m a , as is the case in most of the parameter region giving AN c g < 2 for 
m h ^ 130 GeV. Here we have neglected the masses of the final states and mixing with 
neutral gauginos. Thus, depending on /i and F a , can decay inside the detector while 
leaving displaced vertices. Measuring its decay length would give us information about the 
axion decay constant. 



IV. CONCLUSIONS 



In this paper we have studied the possibility that dark radiation recently suggested by 
observations can be explained by the QCD axion non-thermally produced by the saxion 
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decay. In order to account for AN e s ~ 1, the axion superfield must have a sizable coupling 
to the Higgs sector. We have found that the Higgsino mixing parameter \i is bounded above 
and should be in the range of a few hundred GeV, for the saxion mass lighter than 1 TeV. 
Considering that the saxion mass could be naturally one order of magnitude smaller than 
the soft masses for the MSSM superparticles, the Higgsino can be within the reach of the 
LHC and/or ILC, even if the other SUSY particles are much heavier than 0(1) TeV. This 
will be of great importance especially if the SM-like Higgs boson mass is confirmed to be 
around 124 - 126 GeV (33] . 
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Saxion decay rates and axino production 



In this appendix, we present the partial decay widths of the saxion, and the derivation 
of the axino abundance thermally produced by Higgsino decays. In the decoupling limit 
rriA 3> mw, the saxion decay occurs with 

,2\ 1/2 



¥hh 



mi 1 

l-4-|) A (j, 

m A a 

4 / 2 \ 1/2 / 2 4 

m„ I . mi r \ I . mi, m 



<r-r\ = m —1 1-4— y) ( 1 - 4— - + 12— - ) A a , 

[mi — m h ) A \ m a J \ m A a to* 

22 / 2 \ 3/2 

mimt ( mi N 



[mi - m" h )* \ mi 



r ^// = ^ , 2 ^ - K, (21) 



if the corresponding process is kinematically accessible. Here ky = 2(1) for V = W(Z), and 
Nf = 3(1) for quarks (leptons). The overall factor A CT is defined by 

\B\ 2 \ 2 |/i| 4 



///-, / m% 



A a = AC 2 a B a T a I 1 - ) i^L, (22) 
13 



where we have used the relation m\ = 2\Bfi\/ sin 2/3. On the other hand, if m a > 2/i, the 
saxion decays into a pair of Higgsinos with 

112/ I 1 2 \ 3/2 

r„ M = 8C Xr.|(i-4^) , (23) 

ignoring mixing between the Higgsinos and gauginos. The above decay rate will be reduced 
by the mixing. 

Let us move to the axino abundance. The energy densities of the saxion and thermalized 
radiation are given by 

^)=Po(^) 3 e— >, 

p SM (t) = (1 - B a )T aPo (*X t dt'^e-^'-^, (24) 

where ao = a (to-), po = Po-(to-)- Thermal production of axinos is dominated by the Higgsino 
decays: 

- (cos/^isin/^Jif- (l ' 1 



r, ?2 _ * (eos^sm^§^(l-^) 2 (l + 2 ^), (25) 

for « > and /i ^> m a . The number density of thermally produced axinos is thus 

estimated by 



- = 2J — 5 /_ *a 3 (t)nflo(t) + 



o 



45 / v"^ ^H?^ha,Z 



V (^ <tg + Pt>M) 



+ Z 'Jo dk e^„ + i X M°KW»/)'. 1] j • P6) 

where 2/ = 7//// and z CT = T CT //i. The axino production functions Pj(a;) before and after 
t = t a are approximated by 

f°° ^-13i>2 / ^ \ 9 

e -0.63/x 



/jp /*oo (3 7 2 

d*jf dfc eV ^- A + i «g(x-^ / )x- 2 e- 1 - 26 /', (27) 
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for xq = Max[x, Zf], and 6(x) being the step function. Here we have taken the approximation, 
a(t)/a ~ (t/t CT ) 2/3 ~ (T ff /T) 8 / 3 at t < t a . For z a < 0.2, P t< tA z <r) is always larger than 
Pt>t a { z a)- We also find that the axino abundance computed by the approximated formulae 
is slightly smaller by a factor 2 or 3 than the one obtained by numerically solving the 
Boltzmann equation. However, the above approximation is enough to illustrate how strongly 
the axino abundance depends on T„. 
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